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ABSTRACT 

Recent data from the VVV survey have strengthened evidence for a structural change in the Galactic 
bulge inwards of |/| < 4°. Here we show with an N-body barred galaxy simulation that a boxy bulge 
formed through the bar and buckling instabilities effortlessly matches measured bulge longitude profiles 
for red clump stars. The same simulation snapshot was earlier used to clarify the apparent boxy bulge 
- long bar dichotomy, for the same orientation and scaling. The change in the slope of the model 
longitude profiles in the inner few degrees is caused by a transition from highly elongated to more 
nearly axisymmetric isodensity contours in the inner boxy bulge. This transition is confined to a few 
degrees from the Galactic plane, thus the change of slope is predicted to disappear at higher Galactic 
latitudes. We also show that the nuclear star count map derived from this simulation snapshot displays 
a longitudinal asymmetry similar to that observed in the 2MASS data, but is less flattened to the 
Galactic plane than the 2MASS map. These results support the interpretation that the Galactic 
bulge originated from disk evolution, and question the evidence advanced from star count data for 
the existence of a secondary nuclear bar in the Milky Way. 

Subject headings: Galaxy: structure — Galaxy: bulge — Galaxy: evolution — methods: numerical 



1. INTRODUCTION 

The barred nature of the Milky Way is well-established 
from NIR photometry dS litz & Socrgcl 1991; Dwc k et al.) 
119951: iBinnev et al]|1997t: ,Bissantz fc Gerhard 2002') and 
star counts (iStanek et ai.l 119941: Lopez-Corredo ira et al. 
120051: iBeniami n et al.ll2005l: iCabrera-Lavers et al.ll2007D . 
from comparing HI and CO Iv-diagrams with hydro- 
dynamic models ([Englmaier fc Gerhardl[T999l : lFuxjll99"9l : 
iBissantz et al.l l2003f ), and from dynamical modelling of 
the kinematics of bulge stars (iZhaol 119961 : Wu2 119971 : 
IBissantz et al.ll2004l: iShen et al.|[2oTo[ ). 

One open issue is the structure of the inner Galactic 
bulge, for i?<0.5 kpc. Is this simply the scaled-down, 
higher-density parts of the surrounding boxy bulge? 
Does it become axisymmetric towards the center? Does 
it contain a "classical" bulge component? Or does it even 
contain a secondary nucle ar bar as do about one third of 
barred galaxies like ours ()Erwinll20HT) ? 

This issue is important; for example, a nuclear bar 
misaligned with the main Galactic boxy-bulge bar would 
substantially modify the morphology of the gas flow 
in the central Galaxy, as illustrated by th e simula- 
tions oflRodriguez-Fernandez &: Co mbed (|2008l hereafter 
R FC08) and Namekata et al. (2009). 

'Ni shivama et al.l (|2005 ) analyzed red clump (RC) star 
counts in the inner bulge. Because RC stars are ap- 
proximate standard candles, their magnitude distribu- 
tion i s a proxy for the s tellar density along the line-of- 
sight (IStanek et al.lll99^ . Determining the maximum of 
the RC magnitude distribu tion at v arious longitudes in 
a strip at 6 = — 1°. Nishiva ma et al.l ([2005) found a clear 
change of slope in the RC longitude profile at |/| — 4°, 
separating a steeper slope in the outer bulge from a shal- 
low slope in the nuclear regi o ns. 

Recently, iGonzalez et al.l (|2011aD analyzed new RC 
star counts fro m the VVV survey and found excellent 
agreement with iNishivama et al.l ()2005D . They also de- 



termined the RC maximum versus longitude for b ~ +1° 
and found good agreement with the b = —1° results, as 
expected for a triaxially symmetric structure. This also 
shows that reddening corrections, which differ greatly be- 
tween these latitudes, are unlikely to cause the observed 
change of slope in the longitude profiles. 

The RC star counts therefore indicate a significant and 
fairly sudden change in the structure of the bulge at 
\l\ c=L 4°. Here we analyze the distribution of RC maxima 
for a suitably oriented boxy bulge - barred galaxy model. 
Surprisingly, we find that the RC longitude profiles pre- 
dicted by this model are very similar to the observed 
profiles. We show that this can be traced back to the 
model's density structure becoming more axisymmetric 
towards the center, and finally we compare the asymme- 
try of the predicted star count map in the inner bulge 
region with available analysis of 2MASS data. 

2. A MILKY WAY-LIKE N-BODY MODEL 

The model analyzed in th is paper stems from the 
barred galaxy simulation of iMartinez-Valpuesta et al] 
(|2006( ) and is the same model which we used earlier in or- 
der to argue that the long bar and the boxy bulge in the 
MW are the three-dime nsional and planar components 
of th e MW's main bar (jMartinez-Valpuesta fc GerhardI 
1201 IL hereafter MVGll). It evolved from an initially 
exponential disk with Q = 1.5 embedded in a dark mat- 
ter halo, and developed a prominent boxy bulge through 
a buckling instability after ~ 1.5 Gyr. To resolve the 
nuclear bulge better the particle distribution is sym- 
metrized with respect to the midplane. We consider the 
simulated galaxy at time ~ 1.9 Gyr, some time after the 
instability when the bar has resumed its evolution and 
has regrown through further angular momentum trans- 
fer to the halo. The snapshot chosen for our analysis in 
this paper is that shown in Fig. 1 of MVGll; however, 
we have checked that the results given for the nuclear 
regions below do not significantly change when the snap- 
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Fig. 1. — Distributions of simulated red clump distance moduli 
in different bulge fields (black lines), for an N-body model of the 
Galactic boxy bulge and bar. The model is viewed from a position 
Ro = 8kpc from the Galactic center and Obar = 25° away from 
the bar major axis. The particle positions are convolved with a 
Gaussian luminosity function for the red clump stars with width 
(Tr^C = 0.17 mag. Gaussian fits to the distributions are shown by 
dotted red lines. 

shot m Fig. 2 of MVGll is used instead. 

2.1. Magnitude- dependent simulated star counts in the 
inner bulge 

The steUar population of the Galactic bulge is pre- 
dominantly old 10 Gyr) and has a metallicity 
distributi on with FWHM 1 dex around [Fe/H]~- 
0.1 dex (H rown et all I2O1O0 . In this range the num- 
ber of RC sta rs per unit mass varies within 10% 
(jSalaris fc Gira rdi 2002), so we can assume that the spa- 
tial distribution of RC stars is the same as that of model 
particles (mass). We project the bulge of this model onto 
the sky in longitude-latitude coordinates (Z, h) as seen by 
an observer at i?o = 8 kpc distance from the Galactic 
center (GC hereafter), whose line-of-sight (LOS) is at 
an angle atar = 25° relative to the long axis of the bar 
([Gerhard 2002, MVGll). Dividing the projected particle 
distribution into fields of AZ — A6 = 1°, we sort the par- 
ticles in each field by distance modulus. To account for 
the finite width of the luminosity function for red clump 
giant stars (RC hereafter), we convolve each particle with 
a Gaussian of width ct = 0.17 mag lAlves 2000). The ab- 
solute M/i-(RC) varies by ~ 0.26 mag over the observed 
FWHM metallicity range (.Salaris k, Girardi 2002), cor- 
responding to a (7 of 0.11 mag. Added in quadrature, this 
gives a total a' = 0.21 mag. Using a' instead of tr = 0.17 
mag has a negligible effect on the figures shown below. 

Fig. [1] shows the resulting LOS magnitude distribu- 
tions for fields with 6 = and several I. The main peak 
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Fig. 2. — Maxima of observed and model magnitude distributions 
for red clump (RC) giant stars in bulge fields as a function of 
longitude. Top: Simulated RC maxima for strips with latitudes 
b = [1 ± 0.5°] and [0 ± 0.5°] (black dots), compared with data 
from the VVV survey at b = ±1° IGonzalez et al.ll2011ai , open 
squares). We used Mk = —1.70 to shift the model distance moduli 
to the magnitude scale of the data; see text. Bottom: RC maxima 
in strips with different latitudes, colour-coded as indicated on the 
plot. Individual fields have sizes Al = Ab = 1° and are separated 
in longitude by 1°. The change of slope in the inner few degrees 
seen at low latitudes is absent at 6 = 5°. The horizontal line 
illustrates the assumed distance to the GC. 

corresponds to the inclined bar and boxy bulge; it shifts 
from m-M~ 14.28 at I = 10° on the bar's near side to m- 
M~ 15.14 at I = -10° on the far side. The background 
disk is seen clearly at m-M~ 15.55 in the top panel for 
I = 10°; the slope towards m-M=16 is consistent with 
the disk scale-length. This disk feature is similar in abso- 
lute numbers but weaker in relative numbers for smaller 
\l\ and seems to blend with the far end of the bar at 
I = —10°. A similar phenomenology appears to be seen 
in the VVV data (Gonzalez et al. 2011a, Fig. 2), which 
in addition show a 'background' of K giants not included 
in the Gaussian RC luminosity function. 

We determine the maxima of the simulated RC distri- 
bution in all fields by fitting a Gaussian to the bright- 
est peak of the distribution. Fig. [5] shows the result- 
ing RC longitude profiles for several latitude strips, and 
compares those at 6 = ° and b = 1° to the re- 
sults from IGonzalez et all (|2011aD at b = ±1°. This 
last step requires choosing an absolute K magnitude 
for the RC. To shift the model to the magnitude scale 
of the data, we need to use Mk{KC)= -1.70. The 
most recent Mx(RC)=-1.61 (jAlvcs 200d; iLanev et aH 
I2011f ) for the local RC, co rrected for population ef- 
fects (jSalaris fc Girardill2002D . would predict Mk(RC)~ 
— 1.52. The -1.70 value would thus argue for a shorter 
distance Rq to the GC, wh ile the best cur rent value is 
Ro = 8.3 ± 0.23 kpc (Brunthaler et alllMl . We do not 
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Fig. 3. — Left panel: Face-on surface density of the particles with \z\ < 300 pc, with overplotted maxima of the Hne-of-sight density 
distributions (open circles) and maxima of the simulated line-of-sight RC magnitude distributions (full circles) for particles in the latitude 
range |6j < 2°. Right panel: Same for the particles with 450pc < \z\ < 750 pc, and the latitude range 3° < b < 5°. The long axis of the 
bar is at an angle Obar = 25° relative to the line-of sight to the observer at the assumed galactocentric radius i?o = 8 kpc. Dotted lines 
show directions I = 0, ±2, ±4, ±6, ±8, ±10°, as seen from the observer's position. 



pursue this issue further here. 

Both the observed 6 = ±1° longitude profiles and the 
low-latitude model longitude profiles show the same qual- 
itative result: the steep slope of the RC maxima seen at 
large positive and negative |Z| > 4° flattens for small 
The model profile is slightly flatter than the data inside 
I = ±2° and at ^ > 4°, but steeper around I = -4°. 
The origin of the difference at ? > 6° between model and 
data is unclear; the maxima in the observed distribu- 
tions are very broad and the contribution from the disk 
substantial, but also differential population effects might 
contribute. By contrast with the lower latitudes, the 
6 = 5° model longitude profile shows a markedly weaker 
fiattening in the inner parts |/| < 4°. 

This model was not constructed to match Milky Way 
observations - it is a generic boxy bulge and bar model 
from a simulation; thus we do not expect a perfect match 
with any Milky Way data. The agreement between 
the low-latitude observed and model profiles is there- 
fore even more striking. Previously, the fiattening of 
the observed longitude profiles has been interpreted as 
due to a distinct structure in the inner Gala ctic bulge 
(jNishivama et alj 120051: iGonzalez et al.ll2011a|) . perhaps 
a nuclear bar (see also |Alard| l2001|, RFC08). The simu- 
lated model bulge does not contain such a structure. 

2.2. Origin of the change of slope in the simulated star 

counts 

To understand the cause of the fiattening in the model 
longitude profiles, let us first consider a simple thought 
experiment. We know that the RC maxima for a thin 
bar follow approxima.tely t he major axis of the bar (e.g. 
iCabrera-Lavers et al.ll200"7h . For a thicker bar, e.g., with 
planar axis ratio b/a ~ 0.5, a volume effect arises due 
to the geometry of the cones with fixed Al,Ab. This 
causes the maxima of the RC magnitude distributions 

^ Due to the symmetrization the model star counts at 6 = 1° 
and b = —1° are equivalent. 



to lie behind the maxi ma of the line-of-sight ( LOS) den- 
sity distributions (Cabr era-Lavers et alll2007l MVGll). 
However, as long as the density distribution is approxi- 
mately scale-free, we would expect the RC longitude pro- 
file to have a similar slope as for the thin bar, following 
the orientation of the major axis. 

Now insert an axisymmetric component with high stel- 
lar density in the center of the bar. If its density is suf- 
ficiently high, it will dominate the surrounding bar, and 
since it is axisymmetric, the slope of its longitude profile 
will be nearly zero. By lowering the density of the central 
component, we will therefore be able to arrange a slope 
that is between zero and that of the surrounding barred 
bulge. This simple reasoning shows that the result of the 
previous subsection may be explained if the center of the 
barred bulge has sufficient central concentration and is 
rounder than the outer bar. While if the shape of the 
density distribution were to remain independent of ra- 
dius and only the central concentration were increased, 
in first order only the total number of stars in the central 
LOS would be changed, which does not shift the maxima 
in the histograms. 

Fig. [3] shows the surface density of the barred model 
galaxy in face-on projection, for two slices through the 
model. We see from the left panel that the near-planar 
part of the boxy bulge-bar shows a significant change in 
the axis ratio with major axis length, such that the cen- 
tral parts are nearly round (ellipticity e ^ 0.8) while the 
outer bar is much more elongated (e ~ 0.4). At the same 
time, from the spacing of the contours we see a change of 
slope in the major axis density profile, which sets in at 
approximately major axis length ~ 700 pc corresponding 
in projection to / ~ 4°. The maxima of the LOS density 
distributions and the RC magnitude distributions over- 
plotted on the surface density clearly show the resulting 
change of slope in the inner few degrees. Note also that 
the steep outer slope on the I > side is significantly 
decreased by the volume effect, and that on the I < 
side the outer slope is increased. 
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The Gaussian dispersions found from the longitude 
profile fits in Fig. [T] at I = {-f 0°. -4°, -f °, f °, 4°, f 0°} 
are {0.28, 0.3f, 0.2f, 0.20, 0.24, 0.43} at 6 = 0° (Fig. [U 
and {0.33,0.29,0.21,0.20,0.24,0.40} at & = f°, com- 
pared to the average VVV values at similar longitudes 
and b = ±1° of {0.41,0.33,0.28,0.30,0.33,0.44}. Aver- 
aged over the model longitude profiles for \l\ < 4° and 
b = {0°, 1°}, the average Gaussian (t±4 ~ 0.24 mag, cor- 
responding to a deconvolved FWHM of 0.40 mag or a 
FWHM in the LOS density distribution of ±750 pc. To 
compare with the VVV data, we need to convolve with 
the stel lar p opulation broadening cTpopn — 0.11 mag from 
Section [2. 1[ and with the broadening by residual patchy 
extinction effects within th e 3'-6' subfields used f or de- 
riving the reddening map (jGonzalez et al.ll2011b[ ). An 
upper limit for the latter at 6 = ±1° can be estimated 
by comparing the width in J-K colour of the dereddened 
RC in a field with Ak = 1.5 with a nearby low-extinction 
field at b ~ —1. The observed difference results in 
(Tj_K = 0.28 mag (O. Gonzalez, private communication), 
from which ctk — 0.6ctj_k — 0.17 mag. Adding both val- 
ues in quadrature to the model cr±4 ~ 0.24 mag predicts 
a total a — 0.31 mag, which is close to the observed 
width for the VVV data, a±4 ~ 0.30 mag. The intrin- 
sic RC widths of the model and of the Galactic bulge as 
traced by the VVV RC data may therefore be quite simi- 
lar. For comparison, the projected A/ = ±4° ~ ±560pc. 
Therefore, for both model and data the structure inside 
/ = ±4° is more extended along the LOS than in /. 

By contrast, for the slice with 450pc < |z| < 750 pc 
shown in the right panel of Fig. [3] the change of axis 
ratio in the barred bulge is much less pronounced, and 
also the change of the major axis density slope is less 
pronounced. As a result, the density and RC maxima at 
higher latitudes lie on more nearly straight lines. 

Barred galaxies are often characterized by a rising el- 
lipticity p rofile at consta.nt po sition angle, up to a maxi- 
mum fe.g. lPeletier et al.lll999D, for both disk- and bulge 
dominated galaxies ( Barazza et al.|[2b08[ ). N-body mod- 
els of bars have a range of morphologies; but isodens ities 
have not been studied extensively . Some models in iFuxl 
()1997| ) ; lAthanassoula fc MisiriotisI ()2002| ) appear similar 
to the N-body model used above. 

2.3. Total star count map in the inner bulge 

2MASS star count maps for the Galactic bulge in the 
inner few degrees show a longitudinal asymmetry (Alard 
1200 IL RFC08) which these authors interpreted as inde- 
pendent evidence for the possible existence of a sepa- 
rate nuclear bar. How does the nuclear bulge of our 
simulation compare to these data? Figure |4] shows the 
smoothed model star count map in 13° x 8° and in the 
central 3° x 2°. To construct this map, all particles in 
the simulation have been used, i.e., possible incomplete- 
nesses in the data at the far side of the bulge are not 
modelled. We see that the apparent highest surface den- 
sity is shifted to / ~ —0.3°, as are the centroids of the 
surrounding contours. Also the vertical extent of the iso- 
densities is slightly larger on the I < 0° side for small \l\, 
but becomes larger on the I > 0° side for larger |/|. These 
asymm etries resemble those predicted bv lBlitz fc Spergel 
(|1991f) for the integrated light from the central region of 
an inclined triaxial bulge; however, the star counts here 
are weighted by two powers of distance more than the 
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Fig. 4. — Map of the particle surface density in a longitude- 
latitude plot as seen by an observer 8 kpc from the model galaxy's 
center. Top: central few degrees; bottom: inner bulge region. Con- 
tours are logarithmically spaced. The bar angle «bar = 25° . 

integrated light. 

The pubhshed 2MASS star count maps (|Alardll200l[ 
RFC08) qualitatively show the same asymmetry pattern. 
The most significant difference to the model is that the 
innermost contours in the observed maps are significantly 
more flattened than the corresponding model contours. 
These innermost 2MASS contours are possibly affected 
both by incompleteness (which however would probably 
decrease the apparent flattening) and by blending (which 
might increase it). Similar maps for the VVV data would 
therefore be very valuable. Also, the model used here has 
a somewhat thicker disk than the Milky Way and, most 
importantly, it does not include any dissipative evolution 
such as gas infall and star formation on the orbits around 
the ILR, while s uch star format ion is clearly on-going 
in inner Galaxy (jYusef-Zadeh et al. 2009). Understand- 
ing whether the flattened component seen in the central 
2MASS star count s is consistent with a s ymmetric disk- 
like nuclear bulge ()Launhardt et al.ll2002|) or is related to 
a nuclear bar will require much more detailed data and 
analysis. Here we merely emphasize that the asymme- 
try structure seen in the 2MASS data is not a tell-tale 
evidence for such a secondary bar. 

3. DISCUSSION AND CONCLUSIONS 

We have shown that the RC maximum longitude pro- 
files measured for b = ±1° are well-matched by an N- 
body model of a boxy bulge and bar which formed from a 
bar-unstable disk. This model is the same as in MVGll, 
and we have used the same orientation and scaling to 
avoid parameter fitting. The change in slope of the model 
profiles in the inner few degrees is caused by a transition 
from highly elongated to nearly axisymmetric density 
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contours around this scale. The average FWHM of the 
simulated RC profiles in the range I ~ ±4° (~ ±560 pc 
at i?o = 8 kpc) is ~ 0.56 mag, consistent with the VVV 
data after convolution with estimated stellar population 
and patchy extinction broadening. This corresponds to 
a deconvolved bulge FWHM scale of ±750 pc along the 
LOS in this region. At higher latitudes the transition 
in the model isodensities weakens, predicting that the 
change of slope in the longitude profile disappears. 

The projected star count map of the model has the 
same asymmetry pattern as inferred previously from 
2MASS observations. The central surface density con- 
tours in the Galactic bulge are significantly more flat- 
tened, however. This could imply the existence of a dense 
central disk around the location of the ILR. 



Our results generally support the idea that the main 
part of the Galactic bulge formed from the disk through 
dynamical instabilities. They also show that the evidence 
from star counts so far does not require a secondary nu- 
clear bar in the Milky Way. 
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